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ABSTRACT 

Most direct source detections beyond z ~ 7 are likely to arise from wide-field narrowband surveys of Lya 
emission in the J band. For this to be true, the Lya emission must somehow escape from compact star-forming 
regions (CSR) presumably associated with massive haloes. Since the Lya escape fraction is < 10% from an 
emitting region of size ~ 1 kpc, these objects will be difficult to find and hard to detect, requiring ~ 30 - 100 
hours at each telescope pointing on 8 - 10 m telescopes. For CSR sources, existing large-format IR arrays are 
close to ideal in terms of their noise characteristics for conducting wide-field narrowband surveys where pixel 
sizes are 0.1" or larger. 

However, we stress that Lya can also arise from external large-scale shocks (ELS) due to starburst winds, 
powered by CSRs, ploughing into gas actively accreting onto the dark halo. The winds effectively carry energy 
from the dense, dusty environment of a starburst into lower density regions, where the escape probability for 
Lya photons is greater. ELS emission is expected to be considerably more clumpy (< 100 pc) than CSR 
emission. For ELS sources, IR arrays will need 1-2 orders of magnitude improvement in dark current in order 
to detect dispersed clumpy emission within the environments of massive haloes. These sources require an IR 
camera with small pixels (~ 0.01") and adaptive optics correction {e.g., GSAOI on Gemini), and will therefore 
require a targetted observation of a Dark Age source identified by a wide-field survey. For either targetted or 
wide-field surveys, the deepest observations will be those where the pixel sampling is well matched to the size 
of the emitting regions. We show that there are only 3 -4 J-band windows [z = 7. 72, (8. 22), 8. 79, 9. 37] suitable 
for observing the Dark Ages in Lya; we summarize their cosmological properties in Table 1. 
Subject headings: cosmology; first stars; intergalactic medium; galaxies — high redshift 



1. INTRODUCTION 

In recent years, we have come to recognize that the 'ion- 
ization epo ch' may lie just beyond our current obser vational 
horizon Ccf lFan ef aZ.l2000ll2003HKneib et aZ.l2004l) . Struc- 
ture formation models in CDM cosmolo gies appear to show 
that this epoch took place at z = 7 - 12 (Gnedin & Ostriker 
119971 lHaiman & Loeblfl998l) . Recent wide-angle p olariza- 
tion measurements with WMAP llKogut et aZ.ll20 03) suggest 
the ionization epoch could have been well under way by 
z ~ 17. However, this may conflict with the inferred high 
colum n densities in z ~ 6 quasars (e.g., IWvithe & Loebl 
2004). Since cosmic ionization requires only a tiny fraction 
of the primor dial gas to be converted into stars or black 
holes tLoeb & Barkanall2001[). it is possible to cons t ruct a 
wide range of scenarios JBromrn. Co ppi & Larson J12291 
2002; Nakamura & Umemura 2001; Abel, Brvan & Norman 
2000l I2002t iBromm & Larsod 120041 iBarkana & Loebl 
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Frver. Wooslev & Heger 2001). Clearly, the first observa- 
tions of the Dark Ages will have a dramatic impact on our 
understanding of this new frontier. 

In conventional CDM simulations, CSR sources are associ- 
ated with the cores of massive haloes. Since the distribution 
of dense peaks collapsing o ut of an evolving Gaussian den- 
sity fi eld is well constrained ( Miralda-Escude 2003t IPeacocld 
1999), bright CSR sources are expected to be rare. These will 
be difficult to find even with well optimized wide-field sur- 
veys which exploit large pixels (~0.1"). However, the num- 



ber of detectable sources above a given flux level at a fixed 
epoch is more uncertain. This requires a detailed understand- 
ing of how and when Lya emission is produc ed, and how it 
manages to escape its im mediate environment ( N eufeldll991t 
Haiman & Spaans 1999). Even the most optimistic calcula- 
tions JBaron&Whitei 19871) show that Lya detections will be 
a major observational challenge on 8 - 10m telescopes. The 
possib le discovery of a lensed candidate galaxy at z = 10 by 
Pello et al\ <l2004l) demonstrates the power of lensing to ex- 
tend our observational reach. However, the total number of 
sources accessible this way is limited by the small total vol- 
ume of the universe that is strongly magnified by foreground 
lenses. 

During the Dark Ages, gas accretion onto protogalactic 
cores must have been well under way. Galactic nuclei at the 
highest redshifts observed to date (z ~ 5) exhibit solar metal- 
licities, and theref ore appear to have un dergone many cycles 
of star formation ( Hamann & Ferland 1999). Most galaxy 
cores early in their evolution must have experienced starburst- 
driven winds. Early protogalactic winds will have carried 
large amounts of energy, with relatively low dust/metal con- 
tent, away from the co mplex circumnuclear environment 
(Desiacaue s efa/.l20 04). 

In order to understand physical processes associated with 
the first sources, we will need to resolve the Lya structures. 
But, as we now show, detecting emission powered by outflows 
is a different experimental challenge from that posed by wide- 
field searches. 
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2. THE Lya CHALLENGE 
2.1. Wind-scattered Lya emission 

The difficulty of detecting Lya during the Dark Ages is em- 
phasized in re cent simulations by the GALFORM consortium 
(lLacev et a/J2004tlCole et a/.l2000tl . The GALFORM sim- 
ulations assume that all Lya is produced by star formation, 
and that 10% escapes through the action of galactic winds. 
The simulations are adjusted so as to reproduce the Lyman- 
break and submillimetre galaxy number counts at presently 
observed wavelengths. Down to a limiting flux magnitude 
of f\im = 3 x 10~ 19 erg cm" 2 s" 1 , they predict only about ten 
sources per square arcmin will be observable across the entire 
redshift range corresponding to the J band. The flux limit is 
equivalent to a star formation rate of a few solar masses per 
year which is easily high enou gh to drive large-scale winds 
llHeckmanl20'03llVeilleuxl2003l) . 

The most efficient mechanism for Lya escape is scat- 
tering by neutral H entrained in an outflowing wind 
l lChen & Neufeldlll994l) . The best observed starburst galaxy 
M82 reveals UV-scattered bicones along t he minor axis 
dCourvoisier et a/.lfl990l iBlecha efaZ.lfl990h on a scale of 
500 pc to 1 kpc. This corresponds to a spatial scale of 0. 1 - 

0. 2" at z ~ 8. Thus wide-field searches of rare massive haloes 
will need to exploit pixels on this scale in order to find CSR 
sources. The scattering wind may be clumped but, as we dis- 
cuss below, the individual clumps are likely to fall below the 
detection limit. 

2.2. Wind-driven Lya emission 

We now examine the likelihood that most of the Lya is dis- 
tributed in small clumps which are cooling out of dissipating 
gas, and spread over a larger volume than Lya arising directly 
from star forming regions. 

During the early stages of galaxy formation, collapse is 
likely to be comparable to star formation as a source of energy 
for th e gas. Cold gas accu mulated in the collapse of a small 
halo jRees & Ostrikerl fl977h produces an initial burst of star 
formation. If the energy released into the nascent interstellar 
medium by star formation greatly exceeded that released in 
the collapse, then the bulk of the gas would be unbound from 
the young galaxy. On the other hand, the mechanical energy 
required to limit the initial burst of star formation is compa- 
rable to the binding energy of the gas, which is roughly the 
energy needed to significantly rearrange gas within a galaxy. 
Thus, if the burst of star formation is self-limiting, its feed- 
back cannot fall well short of the binding energy of the gas, 

1. e., the energy released by the collapse. 

Since radiative cooling times of gas at the virial tempera- 
ture are much shorter than dynamical times in low mass pro- 
togalaxies, feedback can only be stored briefly as thermal en- 
ergy. Energy rapidly lost to radiation is ineffective as feed- 
back. Only feedback energy that is converted to kinetic (and 
later potential) energy can be effective in limiting star forma- 
tion on the dynamical timescale. Effective feedback needs to 
induce rapid large-scale flows, i.e., winds, so we assume that 
the primary channel of feedback is through winds. This argu- 
ment requires a moderately high efficiency for the conversion 
of energy released by star formation into wind energy, but 
that is consistent wit h observations of starburst winds {e.g., 
Strickland & StevemfcOOOh . 

If collapse and the initial burst of star formation make com- 
parable energy inputs to the ISM, we can estimate both using 
the spherical collapse model. To form a disk, gas must dissi- 
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FIG. 1 . — Cooling width of a collapse shock at z = 8. This shows the prod- 
uct of postshock cooling time and postshock speed as a function of the shock 
speed, for gas with a preshock electron density of 0.036 cirT 3 and abundances 
of 0, 0.1 and 1 solar (top to bottom). The cooling width is inversely propor- 
tional to preshock density. 



pate at least the vertical component of its velocity dispersion, 
i.e., a 2 /2 per unit mass. For a dissipative collapse — one pro- 
ducing cold gas — the energy dissipated is several times this, 
so we use the estimate of a 2 per unit mass dissipated by gas in 
the collapse. The time taken for a halo to virialize is roughly 
equal to its turn around time, i.e., half of its collapse time, f c . 
If the energy dissipated in the collapse is dissipated in half 
of the virialization time, the dissipation rate 4M g cr 2 /t c , where 
M g is the gas mass involved in the collapse. 

Dark energy is insignificant for collapse at z ~ 8, so that the 
mean density of a halo collapsing at t c is 3ir/(Gt 2 ), giving the 
virial radius for a halo of mass M as R = [GMf c 2 /(47r 2 )] 1 / 3 . 
Treating the halo as an isothermal sphere (GM/R = 2a 2 ), 
then gives a 2 = 0.5(27rGM/f c ) 2 / 3 , so that a ~ 53MJ/ 3 km s" 1 , 
where the halo mass M = 10 10 Mi M and t c = 6.27 x 10 s yr 
[z = 8, in ACDM with (h, O m , ft A ) = (0.7, 0.3, 0.7)]. Little Lya 
can be produced in the collapse unless the gas is heated over 
~ 10 4 K, i.e., unless a > 10 km s" 1 , requiring halo masses ex- 
ceeding 10 8 M Q . 

Using the ba ryon fraction determined by WMAP {fi, = 
ilb/ilm — 0.17; Spergel et a/.l 120031) and assuming that all 
baryons are gaseous in the collapse, the dissipation rate 
in the gas during a collapse is ~ 4f\,Mcr 2 /t c ~ 2 x 

10 40 Mjq 3 erg s -1 . Instantaneous dissipation rates will show 
a significant spread about this value. The fraction emerging 
as Lya photons depends on the shock temperature and the 
fraction of photons that escape. Higher escape fractions are 
favoured by low optical depth and low abundances (low dust 
content). In stand ard ACDM very few 10 10 M halos col- 
lapse before z = 8 JCole et q/l2000l) . 

While the fraction of Lya photons escaping from starbursts 
is small, a significant part of the energy carried off by galactic 
winds is ultimately likely to be dissipated in shocks. These 
are another potential source of Lya photons and, from above, 
the energy in the winds is comparable to that dissipated in the 
collapse. The winds effectively carry energy from the dense, 
dusty environment of a starburst into lower density regions, 
where the escape probability for Lya photons is greater. Star- 
burst wind speeds are high and not strongly dependent on the 
properties of the hosting halo. This means that if starbursts 
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are triggered in halos smaller than 10 8 M©, then the terminat- 
ing shocks of galactic winds may produce Lya although the 
collapse shocks cannot. 

From above, the mean baryon density in a dissipationless 
halo collapsing at z = 8 is 37r/b / (Gmut^) ~ 0.036 cm" 3 , where 
Wh is the mass of hydrogen. Since dissipation significantly 
increases the density of the gas, the typical density of gas run- 
ning into shocks during a dissipative collapse can be some- 
what larger than this. The collapse is fairly chaotic and col- 
lapsing gas is likely to encounter several shocks before finally 
reaching its destination. Provided that the effects of radia- 
tion transfer are not too significant, the depth of the emitting 
region behind a radiative shock is approximately equal to the 
product of the postshock cooling time and velocity. In a strong 
shock, the gas density is increased by a factor of 4, the speed 
is decreased by the same factor and the postshock temperature 
is r ps = 3/imnvl/(16k), where v s is the shock speed and fimu 
is the mean mass per particle. Fig. 1 shows the width of the 
postshock cooling region for a preshock electron density of 
0.036 cm" 3 , for metal abundances of 0, 0.1 and 1 solar. 

Initially, the collapsing gas flow may be coherent on larger 
scales than the postshock cooling length in Fig. 1 . In that case, 
cooling shocked regions are sheetlike and emission is bright- 
est at caustics where we see folds in these sheets projected 
onto the sky. The cooling gas is subject to thermal and other 
instabilities that will generally cause it to fragment on about 
the scale of the cooling length. (The flow into any further 
shocks is likely to be considerably more chaotic.) For haloes 
in the mass range of interest, the shock speeds in the collapse 
are, at most, 100 - 200 km s" 1 . In a pristine gas, the inferred 
size of cooling clumps is roughly 100 pc, or 0.02" at a red- 
shift of z ~ 8. The small size of the cooling region demands 
an infrared imager which utilizes small pixels (~ 0.01") and 
adaptive optics correction. However, we point out that rapid 
winds from starbursts can produce significantly larger cooling 
regions. 

Since gas is compressed significantly in a dissipative col- 
lapse, the time scale for star formation in a clump of col- 
lapsed gas can be much shorter than the dynamical time of 
the collapsing system (~ 10 8 yr). Massive stars take only 
~ 10 yr to produce supernovae, giving plenty of time for 
starbursts and their winds to get underway while other gas 
continues to collapse into the system. Interaction between 
infalling gas and starburst winds can lead to further shock- 
ing of both. This process is observed in M82 where the out- 
flowing starburst-driven wind impinges dir ectly o n infalling 
gas at a radius of 11 kpc iYun. Ho &Lo||l g94|) producing 
observable Ha and x-ray emissio n ([Devine & Ballv Il999t 
iLehnert. Heckman & Weaverll999h . 

3. EXPECTED Lya SURFACE BRIGHTNESS 

We now examine whether wind-powered Lya emission in 
small clumps can be detected with an 8m telescope using an 
adaptive optics imager with small pixels (0.02"). We assume 
that the total energy released in Lya at a fixed star formation 
rate is comparable to what is generated in large-scale shocks 
driven by the starburst, and that this emission escapes without 
attenuation. As discussed in §2.2, this assumes a high con- 
version rate of the su pernova energy into wind energy (e.g., 
Strickland & Stevensl2000l 

From survey s of nearby dwarf galaxies iMartinl 120031 
IVeilleuxl |2()03), we adopt a wind radius of 5 kpc. At 
z ~ 8, the shock occurs at 1" radius and is barely re- 
solved. Naively, if the wind energy escapes along bipo- 
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FIG. 2. — Filter profiles used in our GSAOI calculations. The Lorentzian 
profile (m =1) shown in black is the expected response of a high finesse 
(N=40) tunable filter. The green profile (m — 2) is the low-cost filter option, 
the blue profile (m = 3) is the high-cost filter option. 



lar cones, the Lya surface brightness will be increased 
by the decrease in solid angle compared to a spherical 
wind. This factor can be ~ 10 since most winds appear 
to be highly collimated (jSho pbell & Bland-Hawthorr] 1998: 
Suaai. Davies & Wardl f2003i; However, this factor does not 
apply to a comparison with the expected flux from a CSR, 
since the normal assumption is that the collimated wind ma- 
terial is what renders the central source visible. An escape 
fraction of 10% assumed in the GALFORM simulations is the 
fraction of radiation which escapes along the bicones multi- 
plied by the fraction scattered in the wind. 

For wind-induced emission, in our model, the total power 
available is 10/v lm but this is now spread over many more pix- 
els. The surface brightness of the shock induced Lya emission 
depends critically on how the bipolar shock surfaces are ori- 
ented with respect to the observer. If the bipolar wind lies in 
the plane of the sky, most of the emission is confined to limb- 
brightened arcs over an area of ~ 10 2 pixels depending on the 
curvature of the shock surface. If the bipolar wind is directed 
at the observer, the projected radius of the shocked-induced 
nebula is half the intrinsic radius (60° opening angle). Thus, 
the flux is now dispersed over an area of ~ 10 3 pixels. 

Our expectation is that the Lya flux per pixel will be an or- 
der of magnitude fainter in external large-scale shocks (ELS) 
compared to CSR emission in a large pixel survey (0.1"). We 
believe that it will be important to reach this flux level if we 
are to understand the nature and origin of the Lya emission. 
We have assumed that essentially all of the Lya emission es- 
capes its environs since it is produced in situ. But in the pres- 
ence of any diffuse intergalactic HI component at rest with 
the expanding universe, the shock surface must have a suffi- 
cient redshift in order for the Lya emission to escape to the 
observer. 

We have already noted that CSR emission is rendered visi- 
ble through scattering by neutral H in the wind. From the lo- 



4 



TABLE 1 

Basic cosmological parameters for the four dark J-band 

WINDOWS. 



Wavelength (/ira) 


1.06 


1.12 


1.19 


1.26 


Redshift 


7.72 


8.22 


8.79 


9.37 


Time since Big Bang (Gyr) 


0.66 


0.60 


0.55 


0.51 


Time before present (Gyr) 


12.80 


12.85 


12.91 


12.95 


Physical scale (kpc/arcsec) 


4.92 


4.73 


4.53 


4.35 


Luminosity distance, di (Mpc) 


77283 


82979 


89667 


96398 


Flux, F=L S /Airdl (10" 19 cgs) a 


1.40 


1.21 


1.04 


0.90 


Physical depth, D (Mpc/1000 km s" 1 ) 


1.009 


0.929 


0.849 


0.779 


Comoving depth, D„ (Mpc/1000 km s" 1 ) 


8.798 


8.565 


8.312 


8.078 


Physical volume, 1' X 1' xD (Mpc 3 ) 


0.088 


0.075 


0.063 


0.053 


Comoving volume, 1 ' X 1 ' X D„ (Mpc 3 ) 


58.5 


58.7 


59.0 


59.1 



TABLE 2 

Basic parameters for GSAOI and a DAZLE-like instrument 

OPERATING AT THE GEMINI TELESCOPE. 



Reduced telescope area 
System throughput 
Filter throughput 
Detector pixel size 
Detector read noise 
Detector dark current 
Night-glow continuum 
Source flux 
Source diameter 
Total SNR for source 
Exposure time 



(tt/4)(7902- 1302) 2 cm" 2 

0.27 

0.80 

0.02" -0.1" 
5 e~ pix~' 



2- 5 Rayleigh A 
3 - 30 X 10" 20 
0.04" - 0.20" 

3- 5 
3600 s 



erg cm 



NOTE. — The cosmology is (h,Q m ,Q A ) = (0.7,0.3,0.7). The comoving 
volume is essentially constant over the windows and the total timespan is only 
150 Myr. 

"The expected flux from a source with luminosity Ls = 10 41 erg s -1 , in units 
of 10~ 19 erg cm -2 s _1 , only varies by 50%. 



cal universe, we know that winds arise from highly inhomo- 
geneous starburst cores and entrain most of their mass from 
cooler material in the surrounding disk ( Such kov et all 1996). 
Thus the scattering medium is likely to be highly clumped. 
However, since the total escaping Lya flux in our model is 
fu m , these scattering clumps are likely to be an order of mag- 
nitude fainter than self-illuminating clumps of shocked gas 
induced by the same wind. 

4. SURVEY METHOD 
4.1. Wide-field Lya survey 

We envisage an initial survey with a wide-field, large- 
pixel imager which scans over several fields in order to iden- 
tify an initial list of CSR sources. This is the primary 
motivation of the DAZLE instrument under construction at 
the Institute of Astrono my, University of Ca mbridge for the 
Very Large Telescope (Bland-Hawthorn et al. 2003, see also 
http://www.ast.cam.ac.uk/~optics/dazlei which exploits high 
performance 1Z = 1000 interference filters closely spaced in 
wavelength within the J band. The wide-field (6.9' square) ob- 
servations are built up in interleaved sub-exposures by switch- 
ing between the two filters: the images are then differenced 
in order to detect signals in one band that are not evident 
in the other. This technique has been widely used with the 
Taurus Tunable Filter (TTF) on the AAT and has led to the 
identification of a line-emitting populations out to z ~ 5 (e.g., 
iBarr et a/.l2004T) . 

Each target field with DAZLE will require long expo- 
sures, and may require several fields in order to find a sin- 
gle source. Table 1 summarises the cosmological proper- 
ties of the four 'dark' windows in the J-band origin ally iden- 
tified in the DAZLE design study (Cianci 2003, see also 
http://www.aao.gov.au/dazlei which are clearly evident in 
Fig. 3. 

We now present SNR calculations for a wide-field and an 
adaptive optics near-IR imager on an 8m telescope. The cal- 
culations utilize both airglow and a bsorption spectra at a spec- 
tral r esolution of TZ = 10 000 (cf. lOffer & Bland-Haw thorn 
1998). We adopt instrumental parameters for DAZLE and 
for the Gemini South Adaptive Optics Imager (GSAOI) to 
be commissioned in 2005: the key numbers are listed in Ta- 
ble 2 (McGreg or et a/.l l2003). Note that several of the pa- 
rameters are expressed over a range: the lower value is used 



NOTE. — For parameters shown over a range, the lower bound is specific 
to the ELS (GSAOI) case, the upper bound refers to the CSR (DAZLE) case. 
The one exception is the night sky continuum where we incorporate both 
bounds in both calculations. The reduced telescope area incorporates the 
loss due to the Cassegrain hole (M2 stop). The filter throughput assumes 
an IR detector that cuts off at 2.5/xm; this reduces to 60% for a detector 
cut-off at 5fim due to the need for additional optical density layers for long 
wavelength suppression. The exposure time determines the number of read 
noise contributions to the summed image. 
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FIG. 3. — The calculated total exposure times (in ksec) for the parameters 
listed in Table 2. The solid horizontal lines are the times to reach SNR = 5 
in an initial DAZLE survey. The filled region shows the expected times to 
reach SNR = 3 in a targetted study with GSAOI: the lower envelope is for a 
filter with m = 3, 7^=1000; the upper envelope is for m — 1, 7£=500. The dark 
windows near 1.06/im, 1.19/im and 1.26^tm are the most suitable for Lya 
observations; note how the upper envelope picks these out. The cosmological 
parameters for these four windows are given in Table 1 . 



in our ELS calculation (GSAOI), and the upper value is used 
in our CSR calculation (DAZLE). The GSAOI will utilise a 
HAWAII-2RG array, under development for the James Webb 
Space Telescope, which is expected to have exceptionally low 
dark current. DAZLE exploits a HAWAII-2 array which has 
a much higher dark current, but which is adequate for wide- 
field surveys, as we show below. 

The equivalent luminosity for the quoted Lya flux and the 
physical size of the Lya blobs is given in Table 1 as a function 
of redshift. Since the blob sizes in our calculations are larger 
than the AO-corrected psf of GSAOI, we do not consider the 
J-band Strehl ratio in our calculations. The J-band airglow 
sp ectrum was normalized t o the expected J-band counts given 
in IMcGregor et al\ d2003l) . We consider two limiting cases 
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for the night- glow continuum between the OH lines (e.g., 
lContentl fi996): the actual night-glow continuum level is un- 
known but is bracketed by the quoted values in Table 2. The 
night-glow surface brightness is quoted in Rayleighs (R) per 
Angstrom where 1 R = 10 6 /47r phot cm" 2 s" 1 sr" 1 . 

At a fixed spectral bandpass defined by 7Z, the 
wings of the filter profile must be considered. 
Jones, Bland-Hawthorn & Burton! fl99^> demonstrate 
that the out-of-band blocking of an interference filter with 
m cavities is closely matched to a Butterworth function of 
degree m. We have incorporated the Butterworth profile (see 
Fig. 2) in our calculations. The single cavity (m = 1) is the 
Lorentzian profile of a tunable fi lter in the high finesse limit 
jBland-Hawthorn & Jonesl 119981) . The relatively low-cost 
DAZLE filters utilise m = 2 (two cavity), although there is 
also a high-cost option with m = 3. We assume that placement 
within the instrument does not degrade the effective bandpass, 
an issue discussed at length by Bland-Hawthorn et al. (2001). 

Both m = 2 and m = 3 filters are risky and highly expensive 
items to manufacture since the multi-cavity dielectric coat- 
ings can exceed 10/im in total thickness, requiring hundreds 
or even thousands of layers. As was found in the DAZLE 
study, this can be greatly exacerbated by the need for high 
optical density to achieve long-pass blocking 2 by the interfer- 
ence filter. 

In Fig. 3, the horizontal lines indicate we need 40 ksec per 
field in order to achieve a 5cr detection in one or other band, 
but note that the DAZLE technique surveys twice the volume 
of a single 1Z = 1000 filter image. (If the detection relies on 
the differencing of on-line and off-line bands, its statistical 
significance is reduced to 3.5<7.) An (almost) equivalent strat- 
egy is to halve the exposures per field, and to observe two 
fields in two closely spaced bands. Once tentative sources 
have been identified, in order for the emission to be Lya, the 
source should not be evident in deep ugriz images (i.e., AB 
mag > 28) corresponding to rest-frame Lyman continuum. 

4.2. Targetted Lya survey 

The initial wide-field survey will do little more than iden- 
tify tentative sources for closer study. At this point, we target 
specific sources for further study with a high resolution im- 
ager and adaptive optics. The initial survey will need to target 
fields in the vicinity of IR-bright stars which can be utilised 
for AO correction. A particular advantage of narrowband fil- 
ters is that the AO correction is not hampered by atmospheric 
refraction which can spread point sources over tens of pixels 
within a broadband filter. 

For a targetted study, we require only a single narrowband 
filter, and need only achieve 3er per blob within the Lya neb- 
ula to map the distribution of ionized gas. But Fig. 3 shows 
that the total exposure times per source are 100 ksec. At 
1Z = 1000, there appear to be a dozen useful windows within 
the J- band. However, in our design study for DAZLE ( Ciancil 
2003), we find that only 3-4 windows are practical. For a 
filter placed at the pupil, there is a slight phase effect over the 
field such that the passband at the filter edge is bluer than at 
the centre of the field. For a filter placed in the converging 
beam, the passband is broadened slightly. 

Another concern is the difficulty of manufacturing a 0.1%A 
narrowband filter centred at the prescribed wavelength. Given 

2 An alternative strategy is to exploit an IR array with a sensitivity cut- 
off at a shorter wavelength (A < 2/im). If the dark current can be kept to a 
minimum, we see this approach as preferable to a HAWAII-2RG array. 



the restricted tuning capability of an interference filter, and its 
degraded performance on tilting, it is tempting to consider a 
tunable filter for Lya imaging. However, the upper envelope 
of the filled region in Fig. 3 shows the expected poor perfor- 
mance of a filter profile with Lorentzian wings. 

The above considerations require that we accept a window 
which is actually twice as broad as the design bandwidth. 
In Fig. 3, we note that there is only a handful of good win- 
dows available at 7Z = 500. In terms of photometric stabil- 
ity, the windows near 1.06/im, 1.19/xm and 1.26/im are the 
most ideal. The window at 1.12/xm suffers from variable and 
complex atmospheric absorption. The relevant cosmological 
parameters for these four windows are given in Table 1 . 

The 1.12/im window could be recovered if two IR arrays 
were placed side by side in the image plane, and these devices 
were moved back and forth every few mins, in synchrony with 
nodding by the telescope, such that we build up two separate 
fields of view. This mechanica l form of 'nod & shuffle' (cf. 
Glazebrook& Bland-Hawthorn 2001) would average out at- 
mospheric effects, and at the same time minimize the read 
noise penalty, although it incurs a \/2 penalty from dark cur- 
rent noise. 

If systematics can be controlled, the long total exposure 
times of 30 hr per field are not unreasonable. If there are 
weak OH features in the dark windows, these are likely to 
be variable and may define the systematic noise limit. Note 
that if pixel stability dominates the noise limit, a Dark Age 
experiment cannot be carried out since this leads to short ex- 
posures (Fowler-8 sampling) and results in a huge read noise 
penalty. However, D. B. Hall (2003, personal communica- 
tion) has found that the HAWAII-2RG arrays have compara- 
ble stability to some of the best optical detectors. 

5. THE IMPORTANCE OF DARK CURRENT 

We now show that the IR arrays specified in Table 2 are 
ideally suited to the proposed experiments. If the J-band win- 
dows are truly dark, the read + dark noise contribution defines 
the limiting sensitivity we can reach. This is illustrated in 
Fig. 4 which shows the total exposure time required in the 
1.06/im window to achieve SNR = 3 on a Lya blob with a 
total flux of 3 x 10~ 20 erg cm" 2 s" 1 . Here we assume a fil- 
ter with 1Z = 1000 and m = 3. This calculation is specific to 
GSAOI: the pixels are 0.02" and the read noise is 5 e". At the 
top of the figure, the blobs are spread over many pixels which 
greatly increases the effective background from dark current. 

The results are only weakly sensitive to the number of sep- 
arate exposures since the dark current noise dominates for ex- 
posures longer than about 500 sec. Note in particular how a 
targetted study with GSAOI benefits from a very low dark cur- 
rent array (d.c. < 0.01 e" s" 1 pix" 1 ) for essentially any blob 
size. The HAWAII-2RG array has demonstrably the lowest 
dark current of any existing or planned array technology (D.B. 
Hall 2003, personal communication). 

In Fig. 5, we show the total exposure time required in the 
1.06/im window to achieve SNR = 5 on a Lya blob with a 
total flux of 3 x 10~ 19 erg cm" 2 s" 1 . This calculation is specific 
to the initial DAZLE survey with 1Z = 1000 and m = 3. Note 
how the uncertainty in the dark sky continuum leads to some 
uncertainty in the detectability of Lya blobs. For CSR source 
sizes of 0.2", the existing HAWAII-2 arrays are well suited to 
wide-field Lya surveys of this kind. 

The implication here is that the pixel size must be well 
matched to the expected size of the Lya emitting regions. If 
the pixels are much larger than the emitting blobs, night-glow 
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FIG. 4. — Total exposure times (in ksec) in the 1.06/jm window as a 
function of dark current (horizontal) and blob size (vertical; in arcsec) for 
a detector with 0.02" pixels. We use the parameters in Table 2 specific to 
GSAOI. The exposure times assume SNR=3 per source for a total object flux 
of 3 X 10~ 20 erg cnr 2 s . Note the importance of very low dark current in 
order to detect clumpy Lyct emission. The solid lines assume a high level of 
night-glow continuum between the OH lines (see Table 2), the dashed lines 
assume weak night-glow. 




dork current (e sec pix ) 

FIG. 5. — Total exposure times (in ksec) in the 1.06/^m window as a func- 
tion of dark current (horizontal) and pixel size (vertical; in arcsec). We use 
the parameters in Table 2 specific to a DAZLE-like instrument. The expo- 
sure times assume SNR=5 per source for a total object flux of 3 X 10~' 9 erg 
cnr 2 s . The solid lines assume a high level of night-glow continuum be- 
tween the OH lines (see Table 2), the dashed lines assume weak night-glow. 
For 0.2" source diameters, existing HAWAII-2 arrays are well suited to wide- 
field searches. 



continuum will dominate the exposure and result in exces- 
sively long exposures. If the emitting blobs are much larger 
than the pixels, dark current dominates leading to excessively 
long exposures. 



6. CONCLUDING REMARKS 

It is clear that something quite extraordinary took place over 
a time span of a few hundred million years which led to the 
complete ionization of the intergalactic medium. The epoch 
of (re)ionization is likely to attract a great deal of attention 
over the coming decade. 

We have shown that detecting sources beyond z = 7 is go- 
ing to be very challenging on 8 - 10m telescopes. But since 
30 - 100m telescopes will not be in operation until the next 
decade, narrowband imaging is likely to dominate studies of 
the Dark Ages for years to come. We have outlined a strat- 
egy based on wide-field surveys to identify star formation in 
massive haloes, followed by detailed studies with an adap- 
tive optics imager. The optimal technologies for both kinds of 
study will be available in the next few years. 

It is often observed that major discov eries are mad e within 
five years of a new technology (e.g., iHarwitl [19811) . What 
is less well known is that this statement normally applies to 
experiment s which achieve the systematic noise l imit of the 
apparatus. Glazebrook& Bland-Hawthorn (2001) argue that 
the systematic limit possible with nod & shuffle spectroscopic 
observations argues for total integration times measured in 
months or even years. Since large telescopes are general user 
facilities, it is relatively rare that hundreds of hours are de- 
voted to a single target field. But in order to reach back to the 
Dark Ages, observing programmes of one or two weeks at a 
time will be essential. 
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